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We report the observation of the baryonic B-decay B° — > AtpK~-K + , excluding contributions 
from the decay B° — > At AK~ . Using a data sample of 467 6 BB pairs collected with the BABAR 
detector at the PEP-II storage ring at SLAC, the measured branching fraction is (4.33 ± 0.82 sta t ± 
0.33 syst ±1.13 yl +) x 10~ 5 . In addition we find evidence for the resonant decay B° — > E c (24:55) ++ pK~ 

and determine its branching fraction to be (l.ll±0.30 s tat ±0.09 sys t ±0.29 A + ) x 10~ 5 . The errors are 
statistical, systematic, and due to the uncertainty in the A^ branching fraction. For the resonant 
decay B° -> A+pK*° we obtain an upper limit of 2.42 x 10 5 at 90% confidence level. 
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PACS numbers: 13.25.Hw, 13.60.Rj, 14.20.Lq 

While (6.8 ± 0.6)% [1] 01 a11 -B-mcson decays have 
baryons in their final state, very little is known about the 
decay mechanisms behind these decays and more gener- 
ally about hadron fragmentation into baryons. One way 
to enhance our understanding of baryon production in B 
decays may be to compare decay rates to related exclu- 
sive final states. 

In this paper we present a measurement of the 
Cabibbo-suppressed decay B° — > A+pK~n + 2\. This 
decay can be compared with the Cabibbo-favored de- 
cay B° — > A+pir~ir + , which has been observed by the 
CLEO Q and Belle Q collaborations. The average of 
the branching fraction results from these two experi- 
ments are (12.6 ± 1.3 ± 3.3) x 1(T 4 for B° -> A+ptt-tt+ 
and (2.3 ± 0.3 ± 0.6) x 10~ 4 for the resonant subchannel 
B° — » Z' c (2455) ++ p7r~, where the first uncertainty is the 
combined statistical and systematic error and the second 
one is the error on the Aj — > pK~it + branching fraction. 
If only Cabibbo suppression is taken into account one 
expects the ratio of the corresponding Cabibbo-favored 
and suppressed decays to be close to |Kis/Kid| 2 , where 
V us and Kid are Cabibbo-Kobayashi-Maskawa matrix el- 
ements. Any deviation from this value indicates a con- 
tribution from the additional decay amplitudes possible 
in the Cabibbo-favored decays. 

This analysis is based on a dataset of about 426fb _1 , 
corresponding to 467 x 10 6 BB pairs, collected with the 
BABAR detector at the PEP-II asymmetric-energy e + e~ 
storage ring, which was operated at a center-of-mass en- 
ergy equal to the T(AS) mass (on resonance). In addition, 
a dataset of 44 fb -1 collected approximately 40MeV be- 
low the T(AS) mass (off resonance) is used to study con- 
tinuum background. The BABAR detector is described in 
detail elsewhere For simulated events we use Evt- 
Gen [6j for the event generation and GEANT4 [7j for the 
detector simulation. 

For the decay yl+ — > pK~ir + a vertex fit is performed 
and the invariant mass is required to fall in the interval 
2.277 < m pK7T < 2.295 GeV/c 2 . For the reconstruction of 
the -B-candidate, the mass of the /l+-candidate is con- 
strained to the nominal mass of the At and is com- 
bined with p, K~, and tt + candidates. Afterwards the 
whole decay tree is fitted to a common vertex and the x 2 
probability of this fit is required to exceed 0.2%. 

The selection of proton, kaon, and pion candidates is 
based on measurements of the specific ionization in the 
silicon vertex tracker and the drift chamber, and of the 
Cherenkov radiation in the detector of internally reflected 
Cherenkov light [1| . The proton and antiproton selection 
uses in addition information from the electromagnetic 
calorimeter. The average efficiency for pion identifica- 
tion is about 95% while the typical misidentification rate 
is 10%, depending on the momentum of the particle. The 



efficiency for kaon identification varies between 60% and 
90% while the misidentification rate is smaller than 5%. 
The efficiency for proton and antiproton identification is 
about 90% with a misidentification rate around 2%. 

The separation of signal and background of 
the candidate sample is obtained using two 
kinematic variables, AE = Eg — \fsj2 and 

toes = \J (s/2 + pi • p B ) 2 /E? - \p\ 2 B . Here, y/s is 
the initial center-of-mass energy, E* B the energy of the 
B-candidate in the center-of-mass system, (£"i,pi) is the 
four-momentum vector of the e + e _ system and ps the 
B-candidate momentum vector, both measured in the 
laboratory frame. For true -B-decays toes is centered 
at the B-meson mass and AE is centered at zero. 
Throughout this analysis, _B-candidates are required to 
have an thes value between 5.275 and 5.286 GeV/c 2 . 

After applying all selection criteria there are on aver- 
age 1.16 candidates per event. If the i?-candidates have 
different yl+-candidates we select the one with the invari- 
ant pK~ir + mass closest to the nominal A+ mass If 
the candidates share the same wc retain the one with 
the best vertex fit. 

The significance of the B° — > A+pK~ir + signal is de- 
termined from a fit to the observed AE distribution (see 
Fig. [TJ. As the fit function we use a straight line for 
background and a Gaussian for signal. Fitting between 
-0.12 GeV and 0.12 GeV we obtain 82 ± 17 signal events 
and determine a significance of 8.8 standard deviations 
for this decay. Here, and in the following, we calculate the 
significance as the square root of the difference of 2 times 
the log likelihood of a fit with and without signal compo- 
nent. Like the Cabibbo-favored decay B° — * yl+|>7r~7r + 
the Cabibbo-suppressed decay can proceed via different 
resonant subchannels. Figures and [3] show the side- 
band subtracted vl+7r + and K~n + invariant mass dis- 
tributions, respectively. Here, the signal region corre- 
sponds to \AE\ < 0.024 GeV, and the sideband regions 
to 0.024 < \AE\ < 0.12 GeV. We find evidence for the 
decay B° S c (2A55) ++ pK- (4.3<r) and hints on the de- 
cay B° — > A+pK*° (2.7a). For the determination of the 
significance we use in both cases a second order polyno- 
mial for background and as signal function we use in the 
first case a Gaussian and in the latter a nonrelativistic 
Breit-Wigner. 

For the determination of the efficiency corrected sig- 
nal yield, we divide the phase space into smaller regions. 
In order to account for the resonant substructure, the 
following regions are used: 

1) The i7 c (2455) ++ signal region in the range from 
2.447 to 2.461 GeV/c 2 in m(A+ir+), 

2) the K*° signal region from 0.8 to l.lGeV/c 2 in 
m(K~Tr + ), excluding region 1), and 
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A E [GeV] m (K 7t + ) [GeV/c 2 ] 



FIG. 1: Fitted AE distribution in data with all selection crite- 
ria applied (data points). Shown are all reconstructed events 
B° — > A^pK~-K + . The continuum background, described by 
off-resonance data, is overlaid (histogram). The dashed, ver- 
tical lines indicate the signal- and sideband. 



FIG. 3: Sideband subtracted invariant K n mass with 
the r c (2455) ++ signal region (2.447 < ra(/l+7r + ) < 
2.461 GeV/ c 2 ) excluded. The solid curve is the fit, which is the 
sum of a nonrelativistic Breit-Wigner function and a second 
order polynomial. The dashed curve is the parabolic portion. 
An enhancement at the K*° mass of 896 MeV/c 2 is visible. 
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FIG. 2: Invariant A~£tt + mass in data with the AE sideband 
subtracted. A clear X' c (2455) ++ signal is visible. 



3) all events that are not in region 1) or 2). 

The events in region 3 show no further significant res- 
onant structure, but are also not uniformly distributed 
in phase space. Since we use a phase space model in 
our Monte Carlo simulation we correct the efficiency as 
a function of m(A+pir + ) . We determine the signal yield 
in the different regions by subtracting the extrapolated 
background from the observed number of B-candidates 
in the AE signal region. The background is determined 
with a linear fit to the AE distribution in the AE side- 
bands, 0.024 < \AE\ < 0.12 GeV. For the efficiency es- 
timation we use the same fit strategy as for the signal 
yields, but instead of a straight line we use a second order 
polynomial as fit function to account for the small combi- 
natoric background in the signal Monte Carlo simulation. 



TABLE I: Number of signal events, A s ig, and efficiencies e for 
the three regions used to obtain the signal yield. 



Region 




e 


1 (^c ++ ) 


17.3 ±4.6 


(6.64 ±0.04)% 


2 (iT°) 


26.5 ±9.7 


(8.60 ± 0.07)% 


3 


39.7 ± 12.2 


(8.94 ± 0.25)% 



Here, we use nonresonant Monte Carlo events for regions 
2 and 3 and for region 1 we use B° — * £ c (2Abb) ++ pK~ 
Monte Carlo events since this region is almost saturated 
by resonant events. The number of signal events N S i E , as 
well as the efficiencies e, for the three regions are listed in 
Table |TJ Using these values the overall branching fraction 
is calculated as 



B(B° -► A^pK -7T+) 
1 



(4.33±0.82 stat ±1.13 A +) x 10" 



(1) 



with B(A+ - 
N Bo ± Ar B o_ = 
tion of B°B° 



pK-Ti+) = (5.0 ± 1.3)% [3 and N bE = 
- (467 ± 5) x 10 6 , assuming equal produc- 
and B+B- in the decay of the T(4S). In 
Eq. Jl} and in the following branching fractions, the 
first uncertainty is statistical, while the second one arises 
from the branching fraction of the A£. The final state 
A+pK~Tr + may also include contributions from the de- 
cay B° — > A+AK~. Our cut on the vertex fit probabil- 
ity, however, would strongly suppress this contribution, 
hence the branching fraction ([I]) is understood to not in- 
clude this decay. This is corroborated by the fact that 
the pir + invariant mass distribution shows no A peak. 
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TABLE II: Number of signal events, iV S i g) and the efficiency 
e for the resonant decays via the E c (2455) ++ and the K*°. 

Resonance N s i g e 
r c (2455) ++ f6.0±4.3 (6.15 ±0.04)% 
K*° 20.9 ±7.9 (8.38 ±0.05)% 



For the E c (2455) ++ subchannel we determine the sig- 
nal yield with a ht to the AE sideband subtracted 
m(A+ir + ) distribution. We obtain the signal yield by 
subtracting from the number of events observed in the 
i7 c (2455) ++ signal region the background yield extrap- 
olated from a fit of a second order polynomial to the 
i7 c (2455) ++ mass sidebands. Here, the signal region is 
defined as 2.447 < m(yl+vr+) < 2.461 GeV/c 2 while the 
mass sidebands are 2.426 < m(A+n + ) < 2.447 GeV/c 2 
and 2.461 < m(yl+vr+) < 2.7GeV/c 2 . The efficiency 
is estimated by using the same fit strategy on B° — > 
S c (2455) ++ pK~ Monte Carlo events. Both, the signal 
yield as well as the efficiency for this resonant subchan- 
nel are given in Table [TTJ Using these values we obtain a 
branching fraction of (1.11 ± 0.30 sta t ± 0.29^+) x 1CT 5 
for this subchannel, under the assumption that the 
S c (2455) ++ decays entirely into A+ir + . 

For the K*° subchannel we determine the signal yield 
by a fit to the AE sideband subtracted m(K~ir + ) distri- 
bution, excluding the JC C (2455) ++ signal region. Here, 
we use the sum of a second order polynomial and a 
nonrelativistic Breit-Wigner function in the range from 
0.64 to 1.6GeV/c 2 as the fit function. The nonrela- 
tivistic Breit-Wigner distribution is added in order to 
get a proper background description from the fit. For 
the fit we fix the width and the mean of the Breit- 
Wigner function to its measured values [l|, and deter- 
mine the signal yield by subtracting the integral of the 
background function between 0.8 and 1.0GeV/c 2 from 
the number of events in this region. The efficiency is esti- 
mated applying the same fit procedure to B° — > A+pK*° 
Monte Carlo events. With the obtained values, which 
are listed in Table HH we estimate a branching fraction of 
(1.60±0.61 s tat±0.42 yl +)xl0^ 5 for this subchannel taking 
into account that 2/3 of the K*° decay into K~n + . 

Several sources of systematic uncertainties have been 
investigated. Most of these are derived from studies of 
data control samples and by comparison between data 
and Monte Carlo events. The systematic uncertainties 
arise from the reconstruction of charged tracks (1.4%), 
the charged particle identification (2.4%), and the num- 
ber of BB pairs (1.1%). The uncertainty due to the 
AE background parametrization in data is determined 
by extracting the signal yield with a second order poly- 
nomial instead of a straight line (4.7%). The influ- 
ence of the signal- and sideband definitions is estimated 
by changing their definitions to |A_E| < 0.036 GeV and 
0.036 < | AS | < 0.12 GeV, respectively, and extracting 



the signal yields with these new definitions (3.3%). A 
further systematic uncertainty is the phase space model 
used for the Monte Carlo simulation (1.0%), which is 
determined by reweighting the Monte Carlo events to 
match the observed m(yl+J>7r + ) distribution in data. In 
order to estimate the uncertainties arising from the ap- 
plied m(A+) (3.4%) and \ 2 probability (0.8%) selection 
criteria we vary the criteria by 0.5MeV/c 2 and 0.001, re- 
spectively. The overall systematic uncertainty is 7.5%. 

For the low significance A+pK*° signal, we determine 
an upper limit of 2.42 x 10~ 5 at 90% confidence level. 
This limit is calculated assuming a Gaussian a posteriori 
probability density with a — 0.63 x 10~ 5 which includes 
statistical and systematic errors, and evaluating 90% of 
the integral in the physical region. 

In summary, we observe the decay B° — > A+pK~ir + 
with a significance of 8.8er and measure a branching frac- 
tion of 

B(B° AtpK-TT + ) 

= (4.33 ± 0.82 stat ± 0.33 syst ± 1.13 A +) x KT 5 . (2) 

The ratio of the branching fraction of this decay to that 
of B° -> A+p7r-n+ is 0.038±0.009, which is smaller 
than \V US /V ud \ 2 = 0.0536 ± 0.0020 Q. This is a possi- 
ble indication that additional decay amplitudes for the 
Cabibbo-favored decay are not negligible. Here, and in 
the following the error on the ratio includes statistical 
and systematic uncertainties, while the uncertainty on 
the A+ branching fraction cancels. 

The branching fraction of the decay B° — > 
E c (2455) ++ pK~ is determined to be 

B(B° -> S c (2455) ++ pK-) 

= (1.11 ± 0.30 stat ± 0.09 syst ± 0.29 A +) x 10~ 5 . (3) 

The ratio of this branching fraction to that of B° — > 
£ c (2455)++p7r- HQ is 0.048 ± 0.016, compatible with 
\V US /V ud \ 2 . 

For the decay B° — > A+pK*° the branching fraction is 
determined to be 

B(B° -> AtpK* ) 

= (1.60 ± 0.6W ± 0.12 syst ± 0.42^+) x 10" 5 . (4) 

The 90% confidence level upper limit for this decay is 

B{B° -> A+pK*°) < 2.42 x 10" 5 . (5) 
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